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Abstract 

The Hamiltonian of a charge qubit, which consists of two Josephson junctions is found within 
weh known quantum mechanical procedure. The inductance of the qubit is included from the very 
beginning. It allows a selfconsistent derivation of the current operator in a two state basis. It is 
shown that the current operator has nonzero nondiagonal matrix elements both in the charge and 
the eigenstate basis. It is also shown that the interaction of the qubit with its own LC resonator 
has a noticeable influence on the qubit energies. The influence of the junctions asymmetry and 
the gate capacitance on the matrix elements of the current operator and on the qubit energies are 
calculated. The results obtained in the paper are important for the circuits where two or more 
charge qubits are coupled with the aid of inductive coil. 

PACS numbers: 03.67.Lx, 85.25. Cp, 85.25.Dq, 85.35.Ds 



I. INTRODUCTION 



Josephson-j unction charge qubits are known to be candidates for scalable solid-state 
quantum computing circuits ^,0]. Here we consider a superconducting charge qubit 
which consists of two Josephson junctions embedded in a loop with very small inductance L, 
typically in the pH range. This insures effective decoupling from the environment. However, 
in the practical implementation of qubit circuitry it is important to have the loop inductance 
as much as possible consistent with a proper operation of a qubit. A relative large loop 
inductance facilitates a qubit control biasing schemes and the formation, control and readout 
of two-qubit quantum gates. These considerations stimulated some investigations of the role 
the loop inductance plays in the dynamic properties of charge qubits 0, [5|, j^, where for 
the small loop inductance the corrections to the energy levels due to finite inductance of the 
loop have been found. The corrections have been obtained by perturbation expansion of the 
energy over small parameter j3 = L/Lj, where Lj is the Josephson junction inductance. 

For complex Josephson circuit the construction of quantum Hamiltonian which accounts 
: or finite inductances of superconducting loops can be made with the aid of the graph theory 
3]. This approach has been developed in for systematic derivation of the Hamiltonian of 
superconducting circuits and has been applied for the calculations of the effects of the finite 
loop inductance both for flux Q and charge [l^ qubits. 

In principle, the account for a finite loop inductance (even if it is small) requires for the 
magnetic energy to be included in quantum mechanical Hamiltonian of a qubit from the 
very beginning. It allows one to obtain the effects of the interaction between two-level qubit 
and its own LC circuit. In addition it allows a correct definition of the current operator in 
terms of its matrix elements in a two level basis. 

In this paper we investigate the effect of finite loop inductance and gate capacitance 
for a asymmetric charge qubit, which consists of two Josephson junctions embedded in a 
superconducting loop. 

The construction of exact Lagrangian and Hamiltonian for the charge qubit is given in 
Section II and Section HI, respectively. The approximation for exact Hamiltonian for small 
L is made in Section IV. It is shown that Hamiltonian is decomposed in three parts: qubit 
part, LC-oscillator part and the qubit-LC oscillator interaction part. In this approximation 
the energy levels of the charge qubit are explicitly dependent on the gate capacitance and 
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critical current asymmetry and , in addition, are shifted due to vacuum fluctuations of LC 
oscillator. The current operator both in charge and in eigenstate basis is obtained in Section 
V. It is shown that the asymmetry of critical currents of the Josephson junctions results in 
additional terms in the operator of critical current. The corrections to the qubit energies 
due to its interaction with LC circuit and their dependence on critical current asymmetry 
and on gate capacitance are calculated in Section VI. 

II. LAGRANGIAN FOR THE CHARGE QUBIT 

We consider here a charge qubit in the arrangement, which has been first proposed in ^ 
(see FigHJ. The qubit consists of two Josephson junctions in a loop with very small indue- 




FIG. 1: Charge qubit with inductance coil. 

tance L, typically in the pH range. This insures effective decoupling from the environment. 

As a general case we assume that two junctions have different critical currents Id, Ic2 and 

capacitance Ci, C2. The Josephson energy _Ej = Jc$o/2vr, where $0 = h/2e is the flux 

quantum, is assumed to be much less than the Coulomb energy Ec = {2e)'^/2C, so that the 

charge at the gate is well defined. 

The Lagrangian of this qubit is the difference between the charge energy in the junction 

capacitors and the Josephson plus magnetic energy: 

$2 

L = U - — + Ejicos(fi + Ej2 cos Lp2 (1) 
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where U is the electric energy of JJ's and gate capacities 

TJ = ^ ^ ^ 4- ^ 
2 ^ 2 ^ 2 ' 

$ = J Vidt, where Vl is the voltage drop across the inductance. 
In virtue of Josephson relations 

h 



h 

Ye' 



Vi = —'Pi , ^ = 1,2; V = Vg + —02 



the voltage drop across the inductance is: 



Vl = Vi + V2 



dt dt 



^ {'Pi + 'P2) - 



where $x is the external flux. 

In terms of the phases 992, V^s Lagrangian ^ takes the form: 



(2) 



(3) 



(4) 



L 



+ Eji cos (fi + Ej2 cos (pi 



2(2e)2 



2L 



(5) 



where ipx = 27r$x/'^'o- 

Next we make the known (Likharev and Averin) redefinition of the Josephson phases: 
P^i + 'P2 = P>\ P>\~ P>2 = 25. Lagrangian (jSJ takes the form: 



2(2e)2 V 4 



h h 



— + (Eji + Ej2) cos - cos (5 + {Ej2 - Eji) sin sin 5 
(2e)2 2L ^ ' 2 ^ ' 2 



(6) 



III. CONSTRUCTION OF HAMILTONIAN 



Conjugate variables are defined in a standard way: 

_ 1 _ + C2 + Cg) . - C2 - Cg) . C,K, 

;i9v9 4(2e)2 2(2e)2 ^ 4e 



(7) 



^ 1 ^ ;i(Ci + C2 + c,) • /i(Ci-C2-Cg) . _ cvi^ 

hdb (2e)2 2(2e)2 2e 



(8) 
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From these equations we express phases in terms of conjugate variables: 

Eca 



h 



(9) 



Eca 



(10) 



where Ec = (2e) V2Ce, a = CljC^iC^ + C,), 7 = (C, + - ^2 = ^1 + ^2 + C„ 

Now we construct Hamiltonian: 



H =hn^r,ip + ^1^5(5 — L 



(11) 



Ehminating time derivatives of the phases from with the aid of 0, we obtain the 
final expression for Hamiltonian of the asymmetric charge qubit: 

2 EcOi 



9/ 



2i?j cos — cos b — -Ej^ sin — sin b + -Ej 



2a 



2/3 

where = $o/c/2vr, = (Jci + /c2)/2, ^ = (/c2 - /5 = 27rL/c/$o. 

The first two equations of motion 

IdH IdH 



h dn,, 



are given by Eqs. Q and (fTUI) . Two other equations follows: 

IdH 2Ej if . ^ Ej^ . (f 

ns = -T^TT = ^cos-smd + — 4sm-cosd 

h ob h 2 n 2 



;i2) 



(13) 



n„ 



1 dH 



Ej . f , Ej f . r Ejip-ipx 

■—- sm — cos b H -t cos — sm o 

h 2 2fi 2 hp 



(14) 



Below we consider Hamiltonian (|T2jl as quantum mechanical with commutator relations 
imposed on its variables 



[b, ns] 



(15) 



where = —id /dip, ns = —id/db. 
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IV. APPROXIMATION TO QUANTUM MECHANICAL HAMILTONIAN 



Obviously, Hamiltonian ()12|) is 2D nonlinear oscillator. We assume L is small, so that 
its frequency {LCs) » Ej/h . Therefore we can consider (p as fast variable with fast 
oscillations near the point ipc, the minimum of potential U{ip,6) (see flT^): 

U{ip, 6) = -2Ej cos I cos 5 - Ej^ sin | sin 5 + Ej^-^^^ (16) 

We single out of this potential the fast variable cp, which describes the interaction of the 
qubit with its own LC circuit. 

The point of minimum (pc of U {(p, 6) (fT?)|) with respect to (p is defined from dU /90 = 0: 

V^c = — P sm — cos H cos — sm o (17j 

In what follows we consider 5 as slow variable and expand U{(p, 5) near the point of minimum 
to the third order m. ip {ip = ipc + In the vicinity of ipc the potential U{(p),5) can be 
written as: 

^(^, 5) = Ui^c, ^) + § (l + f ^ ^ + T ¥ ^) ^^^^ 

Ej^3 f . ipC . ^ . r 

(p sm — coso cos — smo 

24^ V 2 2 2 

where (p is the operator conjugate to n^. 

With the aid of (fTTj) we write U{ipc, to the first order in (3: 

U{(pc, 5) = U{ipx, S) = —2E J cos cos 6 — Ej^ sin sin 5 (19) 

Ej(3 f . 2VX 2 . ^ ■ • ox , 2 V^X . 2 X 

' sm cos — - sm (px sm 2o + — cos sm o 



2 V 2 4 4 2 

Therefore, we decompose Hamiltonian ()12|) into oscillator, qubit and interaction parts: H = 
Hose + Hqb + Hint, where 

Hose = Ecanl + - Eca{l - ^)ngn^ (20) 

Hqb = {ns + Ugf - ^^lUg {ns + Ug) +U(v3x,5) (21) 
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Hint = Eca'jn^ns + I cos — cosd + - sm — smdj (22) 

-09 sm cos cos sm o 

24^ V 2 2 2 

r.Ej^2 f . 2 x ^ • • ox , 2V^X . 2x 

+fj—ip Ism -—cos — - sm sm 2o + — cos sm d 
o \ 2 8 8 2 

(sin cos^ 5 — - cos ipx sin 25 — ^ sin (px sin^ 5 



96" V 2 " 4 

In the above equations we disregard the constant term which is proportional to n'^. 

The first term in (j^^ describes the interaction of the phase variables of the qubit, (f and 
6 via the gate, the other terms are responsible for the interaction of the qubit with its own 
LC circuit. 



A. Two- level approximation 



First we quantize (pH) according to {n^ = —ij-', [n^,^] = —i). 



where [a, a+] = 1. 

In addition, we use the two level approximation in the charge basis: ns = | (1 + tz)] 
cos (5 = rx/2; sin 5 = ry/2 with Pauli operators 

r^|0) = -|0);r^|l) = |l); (24) 
rx|0) = |l);rx|l) = |0); 
ty\0) = -z |1) ; ry |1) =z|0). 

In this approximation sin2(5=cos2(5 = 0, since these operators couple charge states which 
differs by two Cooper pairs. Therefore, cos^ (5=sin^ 5=1/2. 

Now we write down Hamiltonian ()2()p2H 0^ within two level subspace in terms of Pauli 
operators tx, ty, tz and oscillator operators a"*", a. 

We obtain the following result: 

H = Wo + Hose + Hgb + Hint (25) 

where 

W^o = ^^^cos^x (26) 
7 



Hose = Ea a^a + - + i 



V2v 



- - {l--f)ng^ (a+ - a) 



64 



sin ifx {a~^ + a) ' 



-?7^ 1 

384^2 V 4 



EcOi 



_ i Eca-f I ^ 



(1 + 2(1 - 7)^3) Tz - TxEj COS ^ - sin ^ 



Tx COS — + Ty- sm — 1 [o^ + a) 



-E 



^ , Tv sm Ty - COS 

96v^ V 2 ^2 



— J (a+ + a) 



where 



Eq 



[2EcEjaV" [2(3EcaV" 



(27) 



(28) 



(29) 



B. The energy levels of the charge qubit 



Here we neglect the interaction of the qubit with its own LC circuit. It is justified if 
13 is sufficiently small so that the energy levels of the qubit oscillator are located much 
higher than the ground level of the qubit. The approximation we make here is to average 
Hamiltonian over the vacuum state, a'^a = 0, of the qubit oscillator. The result is as 
follows: 

1 1 1 

(30) 



where 



H = W + ^Atx + \bty + ]fTz 



W 



32 



Eci- 



2 Ejr,^ 



1 - 



cos (fx 



(31) 



A = -2Ej 



, V \ Vx 

1 cos 

16/ 2 



(32) 



B = -Eji 1 



^ \ . '^x 
— sm 

16 / 2 



(33) 



^ = ^ [1 + 2(1 - l)n,] 



(34) 
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Hamiltonian ()3U|) has the corrections on the order of rf ~ which are due to the vacuum 
fluctuations of the LC oscillator. Since in a charge qubit Ec » Ej these corrections in 
principle might be not very small. 

Hamiltonian (j3(J|) can be made diagonal in the eigenbasis with the aid of the matrix Q]: 

^ /-e-^*sinf cos| \ 

S=\ ' ' 35 

cosf e^^sinf j 

where sin^ = e/AE, cos^ = C/AE, sin^ = B/e, cos^ = A/e- e = VA^ + B^, AE = 

The qubit Hamiltonian in eigenstate basis, therefore, reads: 

S'^HS = W - ^AEaz (36) 

where W is given by (j3T| and 



AE = (l - g)' (cos^f 4=in^f ) + (^) V + 2n,(l-,)f (37) 

As is seen from ()37p the energies of the charge qubit account for its full asymmetry and are 
explicitly dependent on the inductance and on the gate capacitance. 



V. CURRENT OPERATOR 

From the first principles the average current in the loop is equal to the first derivative of 
the eigenenergy relative to the external flux: 

. dE„, 



This expression can be rewritten in terms of exact Hamiltonian of a system: 



(38) 



, , dH , , 

/ = (-1 ^ \n) (39) 
From ()39p we would make ansatz that the current operator is as follows: 

I = 1^ (40) 

However ()4()|1 is not a consequence of ()39|1 . Therefore, the ansatz ()4()|1 must be proved in 
every case, since the current operator in the form of Eq. ()40|) has to be consistent with its 
definition in terms of variables of Hamiltonian H. The prove for our case is given below. 



The current operator across every junction is a sum of a supercurrent and a current 
through the capacitor: 

ii = Ici sin if^ + ^Ci<fi (i = 1, 2) (41) 

We are interested in the current through the inductance coil, Ji in ()4H1 (see FigGJ). 
Direct calculation of Ji with the aid of ()9pi()|14|13|l yield the result: 



c- 



(42) 



which is independent of parameters of a particular junction in the loop. From the other 
hand the expression can be obtained from our Hamiltonian (fT^ with the aid of ()4Up . 
Therefore, the equation (^Ul) gives us the true expression for the current operator. It is 
important to note that the proper expression for the current operator (jl^ cannot be obtained 
without magnetic energy term in the original Lagrangian (pj. 



7^ 0. Therefore, an eigenstate of H cannot 



It follows from ^ and (021) that 
possess a definite current value. 

For the charge qubit the current operator can be obtained from ()30p with the aid of its 
definition (gU)): 



2n dW 



7T 



--TX + -Ary 



$0 $0 

The transformation of in the eigenstate basis yeild the result: 



(43) 



S ^IS = Io + Izoz + ixox + iyoy 



(44) 



where 



27r dW 
$0 (9$x 



(45) 



2 9$ 



X 



TT 



E 



J 



$0 AE 



e 

4 



16 



sm (fx 



(46) 



-Ej 

$0 



1 _ !L 

16 



sm 



cos2 ^ + ^ sin2 ^ 



— Ej- 
$0 2 



1]^ 
16 



cos ■ 



^x 



1 + ^ ^- 



1 ) sin^ £f 



COS 



2 



+ ^ sin2 £f 



(47) 



(48) 
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Therefore, the current operator is not diagonal neither in the charge basis no in eigenstate 
basis. If we neglect the inductance and the current asymmetry {P = 0, = 0), we obtain: 
/y = 0, 

^ - sint^x (49) 



<l>o AE 

TT C . Cpx 

J-x = - 

where 



'x = -—Ej—- sin — (50) 



AE = ^AE] cos2 ^ + C2 (51) 

The existence of nondiagonal elements of the current operator in eigenstate basis is important 
if we consider the inductive coupling of several qubits. 



VI. INTERACTION OF THE CHARGE QUBIT WITH ITS OWN LC CIRCUIT. 
CORRECTIONS TO THE QUBIT ENERGIES 

Here we enlarge the Hilbert space to add to two qubit states two photon states of LC 
resonator, a^a = 0,1. The transformed Hamiltonian ()25|). which accounts for transitions 
between ground and excited state of LC resonator will read: 



where 



S-^HS = W + Pa+a + Qi (a+ ~ a) + Q2 (a+ + a) 



--AEaz + R [a^ 



i) (Caz + Aax - Ba^ 



+S (a+ + a) {Zaz + Xax + Yay) + Ta+a {^z 



AC 



BC 



-Ox 



32 



^2 / ^2 
1 + — - I 1 - — I cos v^x 



(52) 



(53) 



g.^.^[7-2(l-7) 



(54) 



Q2 



128^2 



sm (px 



(55) 



R 



i Eca'j 
2y/2 r]AE 



(56) 
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S = (57) 



16 



AE \ A 



X = B\ l + ^ ^ I (59) 

Y = a\i + ^ ^-^ I (60) 

^ = 7^ .\ (61) 

The operators ax^cry^crz are defined on the eigenstates | t), | i): 

crz|T) = IT); crz|i) = -|i); 

^x|T) = |i); ^x|i) = IT); (62) 
^y|T) = '^li); <yY\i) = -^\]) 

In addition, we restrict photon subspace to photon numbers n=0,l. The basis set for our 
photon+ qubit system is: 

|OT) = |0)®|T); |Oi) = |0)®|i); ^^^^ 
|1T) = |1)®|T); |li) = |l)®|i) 

Within this basis the wave function for Hamihonian ()52|1 is decomposed as: 

^ = a|0 T) + 6|0 i) + c|l T) + i) (64) 
The Schrodinger equation = E^^ takes the form: 

[a {W - \AE - E) +c{Q2-Qi- RC + SZ) -d{RA + iRB - SX + iSY)] |0 t) 

+ [b{W + \AE -E)+c {-RA + iRB + SX + iSY) + d{Q2 - Qi + RC - SZ)] |0 i) 

+ [a {Qi + Q2 + RC + SZ) + b {RA + iRB + SX - iSY) 

+c{W + P- \AE + T-E)~d^{A + iB)] |1 t) 

+ [a {RA - iRB + SX + tSY) + b{Qi + Q2- RC - SZ) 



-c-^ {A-tB)+d{W + P + \AE - T - E)] |1 i) = 



£2 



(65) 
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The energy levels are defined by equating of the determinant of this equation to zero. In 
order to simplify the problem we assume the inductance of the qubit is very small. In this 
limit we may put W = 0, P = Eq, Q2 = 0, S = 0, T = 0. The Hamiltonian is reduced 
to: 

H = Eoa+a + Qi (a+ - a) - ^AEaz + (a+ - a) {Caz + Aax - Bay) (66) 

It is important that the inductance cannot be eliminated at all in the two photon approx- 
imation, since the quantity f3 is in the denominators of Eq and R. The Schredinger equation 
for Hamiltonian (jMj) is as follows: 

[a (-|AE -E)+c {-Qi - RC) ~dR{A + iB)] |0 t) 
+ [b {\AE -E)+cR i-A + iB) + d {-Qi + RC)] |0 i) 
+ [[a{Qi + RC) + hR{A + iB) + c{E^~\AE-E)] |1 t) 
+ [aR {A -iB) + h {Qi ~RC) + d {Eq + \AE - E)] |1 i) = 

The energy levels are defined by equating of the determinant of the following matrix to zero: 

-\AE-E -Qi~RC -R{A + iB) 

\AE-E -R{A-iB) -Qi + RC 

Qi + RC R{A + iB) Eo-^AE-E 
R{A-iB) Qi-RC Eq + \AE-E 

Below we calculate the energies for the following set of the qubit parameters: Ej = 
4.64 X 10~^^J, Ec = lO-Ej, Ug = —0.5; ^ = 0.1. The calculations have been performed for 
two cases. For noninteracting qubit we used the expression (jHTj) . where the LC circuit only 
renormalizes the energy due to the vacuum fluctuations of LC ocsillator (factor rj in (jSZ|))- 
For the qubit which interacts with its own LC circuit we solved the equation (j68|) . 

The plots of the qubit energy levels are shown on FigI21 As is seen from these plots the 
finite value of f3, though it is rather small, modifies the gap between ground and first excited 
energy levels. The calculations show that at the point (px = ^r, where the gap is minimum, 
the gap for noninteracting qubit is AE = 0.14ii^j. The interaction of the qubit with its LC 
circuit reduces the gap to AE = 0.11-Ej. 

The dependance of minimum gap (at the point ng = —0.5, (px = vr) on /5 is shown on 
FigEl As is seen from FigOthe noninteracting qubit is slightly modified by the inductance. 
A small reduction of the gap with the increase of (3 is due to the factor rj in (jHTj) . However, 



= (68) 
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■2,0 



T 




T ' r 

2 4 



T 
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FIG. 2: The plots of qubit energies vs magnetic flux <j)x for 7 = 0.01, /? = 0.001. The red curves 
(a) are the ground and excited states for a qubit which does not interact with its LC resonator. 
The blue curves (b) are the same as (a) but for a qubit which interacts with its LC resonator. 
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FIG. 3: The dependance of minimum gap on (3. The black boxes and red circles are calculated for 
noninteracting qubit. The blue triangles are those for interacting qubit. 

the reduction of the gap with the increase of (3 for interacting qubit is much more significant 
(black triangles on Figl^l). It is important to note that this effect is more pronounced for 
relative large 7's. For 7 = the interaction with LC circuit does not alter the energies 
(compared to those for noninteracting case). 
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